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A PRELIMINARY INVESTIGATION OF EMBRYONIC 
INBREEDING DEPRESSION IN TWELVE 
SPECIES OF LEPIDOPTERA 
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ABSTRACT. Inbred (sib-sib) matings were made in population cultures of twelve 
species of Lepidoptera ( Boloria bellona, B. selene , Phyciodes tharos , P. batesii , P. 
campestris, Junonia coenia , Pararge aegeria, P. megera , Pieris rapae, Arctia caja, and 
Cisseps fulvicollis) and the resulting loss of egg hatchability was compared with con- 
trol, outbred matings in the same cultures. A relatively extensive series of inbred mat- 
ings within a single population of Phyciodes tharos revealed that embryonic genetic 
load is composed of very few, strongly deleterious genetic units. This seems to be true 
of all of the lepidopteran species tested. Genetic load values for nine of the species 
seemed homogeneous and varied between 0.5 and 2.0 lethal equivalents per zygote. 
There was only a single value for the tenth species; this may be anomalous. Two 
species of the group seem to have significantly lower and higher genetic loads than 
the other species. This may be a result of differences in population structure. 



Genetic variation within populations and species theoretically has 
a highly adaptive role in enabling the organism to respond to varying 
environmental stress (e.g., Dobzhansky, 1970). However, a portion of 
this variation (the “genetic load”) consists of genes that, when ex- 
pressed, result in a loss of fitness for the carrier. Since segregation 
and recombination in sexual, non-selfing organisms lead inevitably to 
the production of disadvantageous homozygotes, each population pre- 
sumably must strike a balance between population size and structure 
(and thus levels of inbreeding) on the one hand and size of the genetic 
load on the other. Under conditions such as the laboratory, where 
close inbreeding suddenly becomes the rule, genetic load is ex- 
pressed as inbreeding depression. Comparative studies on inbreeding 
depression in different species and populations within species should 
provide information on the composition of genetic loads and even- 
tually enable some general inferences to be drawn on the species’ 
natural levels of inbreeding and population structure. This paper rep- 
resents a preliminary attempt to estimate genetic loads for a small, 
heterogeneous group of lepidopteran species. 

The magnitude of the genetic load in a population is usually esti- 
mated by observation of the reduction from normal survivorship dur- 
ing a given sensitive period (i.e., during a time of rapid tissue devel- 
opment and growth) in development of progeny from a series of 
matings of a known degree of inbreeding. This approach has been 
used for populations of Drosophila (e.g., Dobzhansky et al., 1963), 
Douglas-fir (Sorenson, 1969), humans (Morton et al., 1956), and do- 
mestic animals (Pisani & Kerr, 1961; Sittmann et al., 1966). A second 
method has been used for Drosophila. This uses genetic markers as 
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a means of estimating the frequencies in a population of chromosomes 
bearing lethal and sublethal genes. When both methods are used si- 
multaneously on the same Drosophila populations (e.g., Dobzhansky 
et ah, 1963; Malogolowkin-Cohen et ah, 1964) the former method 
gives a consistently lower result, indicating synergism among the 
components of the load. Because of this synergism, estimates of load 
magnitude made under conditions of intense inbreeding can be con- 
sidered only as relative values. However, the intensity of inbreeding 
in the present experiments was less than that necessary to produce 
significant synergism (Kosuda, 1972) and should reflect actual levels 
of genetic load. 

During the course of breeding work involving maintenance of 
twelve species of Lepidoptera ( Boloria bellona bellona Fab., B. se- 
lene Schiffermiiller (ssp. mijrina Cramer and sabulocollis Kohler), 
Phyciodes tharos tharos Drury, P. tharos pascoensis Wright (proba- 
bly best considered a separate species (Oliver, in prep.)), P. batesii 
Reakirt, P. campestris Behr (ssp. montana Behr), Junonia coenia 
coenia Hiibner (all Nymphalidae); Pararge aegeria L. (ssp. tircis But- 
ler), P. megera megera L. (Satyridae); Pieris rapae rapae L. (Pieridae); 
Arctia caja caja L. (Arctiidae); and Cisseps fulvicollis Hiibner (Cte- 
nuchidae)) at various times during a period of ten years, I took the 
opportunity of making observations on control and inbred broods 
reared under the same conditions and usually simultaneously. I have 
made estimates of the magnitude of genetic load in each population 
culture from a comparison of survivorship of progeny from control 
and inbred matings. In these species by far the greatest expression of 
genetic load occurs during embryonic development (Oliver, unpub. 
data). Data on embryonic survivorship is also much easier to gather 
than that on later life stages. For these reasons this paper considers 
only embryonic viability among the twelve species. 

Unfortunately, there is little hard data on the population structure 
of these species. P. t. tharos has a wide range of habitats and is gen- 
erally abundant; there is good evidence for near panmixis over the 
eastern United States (Vawter & Brussard, 1975). P. t. pascoensis is 
more habitat-restricted and tends to show significant population dif- 
ferentiation (Oliver, in prep.); P. c. montana probably has a similar 
population structure. P. batesii in the northeastern United States oc- 
curs in small, highly isolated populations (Oliver, 1979). My field 
observations indicate that P. megera (Oliver, 1972b) and P. aegeria 
show little individual motility and tend to occur in localized popu- 
lations, as do B. selene , B. bellona, and C. fulvicollis. A. caja is com- 
mon and general in central Europe; P. rapae is abundant in almost 
all open areas in the eastern United States. J. coenia is a highly vagile 
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species that tends to form large colonies which may often be founded 
by single fecundated females (e.g., Shapiro, 1978). 

Almost nothing is known about the composition of the genetic load 
in Lepidoptera. If the heterozygosity involved in genetic loads is sim- 
ilar to the great genic variation that has been observed in natural 
populations of most organisms (Lewontin, 1974; Richmond, 1972), 
including butterflies (Burns & Johnson, 1967), then the chances are 
very small of picking at random even two parents having a level of 
variation differing significantly from the population average. The data 
of Dobzhansky et al. (1963) seem to support this view; no correlation 
was found between viability and pedigree in a long series of inbred 
Drosophila cultures. However, when lethals or semilethals are in- 
volved, very high levels of inbreeding depression may be due to the 
actions of only a few genes. 

Methods 

One to six cultures were maintained for each species. Detailed in- 
formation on culturing techniques has been given in earlier papers 
(Oliver, 1972a, 1972b, 1977, 1978, 1979). Each culture was begun 
from parents caught in a separate locality. Embryonic survivorship in 
the progeny from inbred matings between siblings (coefficient of in- 
breeding, F, = .250) was compared with that from control, outbred 
intrapopulation matings (“combined controls” in Table 1). In each 
case the siblings used were the progeny of an outbred (usually wild) 
intrapopulation mating having normal hatchability. 

Each of the cultures was begun with a comparatively small number 
of wild parents, from 2 to 22. Because of uncertainty regarding the 
reliability of data from very small population samples, the program 
was carried out in two parts. The first (Part A) consisted of observa- 
tions on cultures of all species maintained from 1968 to 1978, the 
second (Part B) of a more systematic set of inbred matings made dur- 
ing 1973 and 1975 using a population of Phyciodes tharos from west- 
ern Pennsylvania. In this part six isofemale families were established 
and 9 to 14 successful sibling pairings made within each family. An 
additional two families (73-23 and 73-31) were derived by outcrossing 
the Fi progeny of the original wild parents of these families. Sib mat- 
ings were then made within each of these two additional families. 

Cultures were begun from individuals caught at the following lo- 
calities: 

Boloria b. bellona — MASSACHUSETTS: Acton, Middlesex Co.; Pennsylvania: Forward 
Twp., Allegheny Co. 

B. selene myrina — Massachusetts: Acton, Middlesex Co. 

B. s. sabulocollis — South Dakota: Custer Park, Custer Co. 



Table 1 . Mean hatchabilitv of eggs and embryonic genetic loads (L. E./zygote) in control and inbred lines of Lepidoptera. Coefficient 
of inbreeding (F) = 0 for controls (c); (F) = .250 for inbred (i) broods. 
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N.W. Germany (c) 5 ? 4628 .936 ± .055 .888 ± .104 

N.W. Germany (i) 11 ? 10,587 .861 ± .159 .522 ± .287 2.125 

1 Showing visible evidence of embryonic development. 

2 Mean ± standard deviation given if no. of broods >1. 



56 



Journal of the Lepidopterists’ Society 



Phyciodes t. tharos — Florida: 4 mi. E of Cedar Key, Levy Co.; Texas: San Antonio, 
Bexar Co.; West Virginia: Spruce Knob, elev 4500 ft, Pendleton Co.; Pennsyl- 
vania: Upper Tyrone Twp., Fayette Co. 

P. t. pascoensis — Montana: 1 mi. S of Hall, Powell Co.; Alberta, Canada: 6 mi. E 
of Nordegg, Red Deer. 

P. campestris montana — California: Lang Crossing, elev 4500 ft, S. Fork Yuba R., 
Nevada Co. 

P. batesii — New York: Syracuse, Onondaga Co. 

Junonia c. coenia — Florida: Wauchula, Hardee Co. 

Pararge aegeria tircis — England: Pamber Forest, Hampshire. 

P. m. megera — England: Oxford, Oxfordshire; France: Boulogne-sur-Mer. 

Pieris r. rapae — PENNSYLVANIA: Upper Tyrone Twp., Fayette Co. 

Arctia c. caja — Northwest Germany. 

Cisseps fulvicollis — CONNECTICUT: Woodbridge, New Haven Co.; MASSACHUSETTS: 
Littleton, Middlesex Co.; TEXAS: San Antonio, Bexar Co.; IOWA: Dubuque, Du- 
buque Co.; South Dakota: Custer Park, Custer Co. 

Observations on embryonic survivorship included counts of eggs 
showing evidence of embryonic tissue differentiation (a darkening of 
egg color from pale yellow or green to brown or black) and of eggs 
hatching. An increase in failure to show embryonic development in 
the inbred broods was taken as evidence of very early embryonic 
mortality. There is a possibility that there was lowered fertilization 
rate of eggs in the inbred matings, but the males used in these matings 
were of the same age and condition as those used for laboratory con- 
trol matings. Eggs fiom these control matings showed fertility equal 
to that from wild-fecundated females. Both visible fertility and inci- 
dence of egg hatch are combined under the term “hatchability.” 

The Wilcoxon Two-Sample Statistic (Owen, 1962) was used to com- 
pare differences in hatchability of broods in Part B. Calculations of 
load magnitude were made using the equations presented by Freire- 
Maia (1964), which are a simplification of those of Morton et al. (1956). 
Since the comparison is between inbred and control broods with F = 
0, it is possible to use Freire-Maia’s equation 

B = log(Sj/S c ) 

-0.4343Fi 

where B is the magnitude of the genetic load in lethal equivalents 
per zygote; S* and S c denote the mean survivorship of the inbred and 
control broods, respectively; and Fj is the coefficient of inbreeding of 
the inbred group. Morton et al. (1956), who used the term “mutant” 
to refer to any gene deleterious when homozygous, defined lethal 
equivalent (L. E.) as “a group of mutant genes of such number that 
il dispersed in different individuals, they would cause on the average 
one death, e.g., one lethal mutant, or two mutants each with a 50 per 
cent probability of causing death, etc.” 
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Results 

The results of Part A are shown in Table 1. In every species in- 
breeding caused a significant reduction in hatchability, but the 
amount varied greatly from one population and species to another. In 
the control broods hatching failure was due in the great majority of 
cases to infertility rather than to embryonic inviability. In the inbred 
broods, however, there were usually substantial reductions in both 
the incidence of visible embryonic development and in late embry- 
onic viability. There was interspecific and interpopulation variation 
in the inbred broods in the proportion of mortality occurring before 
and after visible evidence of embryonic development. Variation in 
inviability was many times higher for the inbred broods than for the 
controls. The amount of variation was not a function of the magnitude 
of inbreeding depression. There were large differences between 
species and between populations in the size of the genetic load (L. 
E. /zygote calculated from total hatchability) and of the standard de- 
viation in viability observed for each series of broods. 

The results of Part B are shown in Table 2. Hatchability within 
each family line of this population of Phyciodes tharos varied greatly, 
but only the family showing the least inbreeding depression (75-3) 
differed significantly (P < .005) from those showing the greatest (73- 
23, 73-31, 75-4). 



Discussion 

The degree of interbrood variation in hatchability (indicated by 
standard deviation) shown within each family line of P. tharos and 
the other Lepidoptera used here reveals that genetic loads in these 
species consist mainly of relatively few genetic units with strongly 
deleterious effects when homozygous. These genetic loads do not 
seem directly to reflect the genic variation that has been shown to 
occur in Lepidoptera. 

Fig. 1 shows the distribution of L. E. /zygote values for the 85 inbred 
P. tharos broods summarized in Table 2. The graph consists of four 
main peaks, with values clustered at or near 0, 1, 2, 3, and 4 L. E./ 
zygote. A distribution of this type could result from the action of only 
several strongly deleterious (i.e., lethal) recessive genes in each fam- 
ily, perhaps along with a few much weaker sublethals. 

In Table 2 there is a considerable spread in L. E. values, even 
though only the very lowest value is significantly different from the 
highest. Thus, there is a spread on either side of the mean of about 
0.7 L. E. /zygote that is not significant here. If we assume for the 
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Table 2. Mean embryonic genetic loads expressed in Fj progeny of sib-sib matings 
in 8 lines of Phyciodes tharos from Upper Tyrone Township, Pennsylvania. 



Family line 


No. of 
broods 


No. of 
wild parents 


No. of 
eggs 


Hatched/laid 1 


L. E./zygote 


controls 


12 


24 


3037 


.993 ± .011 




73-1 


14 


2 


5884 


.765 ± .154 


1.043 


73-4 


9 


2 


2459 


.803 ± .184 


0.850 


73-23 


9 


4 


2663 


.596 ± .272 


2.041 


73-31 


9 


4 


2595 


.591 ± .146 


2.076 


75-1 


12 


2 


3749 


.779 ± .132 


0.971 


75-2 


11 


2 


3046 


.683 ± .172 


1.494 


75-3 


12 


2 


3021 


.888 ± .113 


0.448 


75-4 


9 


2 


2571 


.671 ± .190 


1.623 
Jc = 1.318 



1 Mean ± standard deviation. 



moment that this amount of variation applies also to the other species 
in Table 1, all except P. batesii ,J. coenia, and C . fulvicollis probably 
fall within the normal range of P. tharos and form a fairly homoge- 
neous group. The single value for P. batesii falls at the upper limit 
of values for P. tharos (above 98.81% of values in Fig. 1), but this 
single sample cannot be considered significant. All three values of J. 
coenia are at the lower level for P. tharos; there is a less than .001 
chance of any three randomly picked values of P. tharos from Fig. 1 
all falling this far to one end of the graph. The overall mean L. E ./ 
zygote value for C. fulvicollis falls at the upper end of the graph. 
From this it seems likely that C. fulvicollis and at least this popu- 
lation of/, coenia carry genetic loads that are, respectively, larger and 
smaller than those of the other species. A very small genetic load 
would seem adaptive for the inbred populations resulting from the 
sort of colonization events that occur in J. coenia. Too little is known 
about the population structure of C. fulvicollis to speculate on the 
comparatively large genetic load it carries. 

The lethal equivalent values calculated here are generally in line 
with those estimated for other insects. One to two L. E./zygote have 
been estimated (by survival from egg to adult) for Drosophila (Dip- 
tera) (Dobzhansky et al., 1963; Stone et al., 1963; Malogolowkin-Co- 
hen et ah, 1964) and for Tribolium (Coleoptera) (Levene et ah, 1965). 
for humans a value of 3 to 5 L. E./zygote (by survival from birth to 
sexual maturity) has been established by Morton et ah (1956), whereas 
embryonic loads in Douglas-fir average about 10 L. E./zygote (Sor- 
enson, 1969). In none of these organisms has a definite relationship 
been shown between magnitude of load and population structure. 

1 he general similarity of load magnitude among the very different 
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Fig. 1. Distribution of embryonic genetic loads in lethal equivalents per zygote in 
85 inbred broods of Phyciodes tharos from Upper Tyrone Township, Pennsylvania. 



insects thus far examined may indicate that genetic loads are more 
related to some internal genetic balancing responsible for the meta- 
bolic integration of the individual organism than to population struc- 
ture. However, there are surely situations (as in species colonizing by 
single fecundated females) where adjustment of load magnitude is 
highly adaptive. 



Acknowledgments 

I am grateful to C. L. Remington, Thomas Uzzell, and Wyatt W. 
Anderson for discussions of the early stages of this work, and to E. G. 
Leigh, Bruce Wallace, and the late Theodosius Dobzhansky for their 
very helpful suggestions and comments on various versions of the 
manuscript. Portions of this research were supported by funds pro- 
vided by the Department of Biology, Yale University, and by a Grant- 
in-Aid of Research from the Society of the Sigma Xi. 

Literature Cited 

Burns, J. M. &: F. M. Johnson. 1967. Esterase polymorphism in natural populations 
of a sulfur butterfly, Colias eurytheme. Science 156: 93-96. 

DOBZHANSKY, Th. 1970. Genetics of the Evolutionary Process. Columbia Univ. Press, 
New York. 505 pp. 

, B. Spassky, & T. Tidwell. 1963. Genetics of natural populations. XXXII. 

Inbreeding and the mutational and balanced genetic loads in natural populations 
of Drosophila pseudoobscura. Genetics 48: 361-373. 

Freire-Maia, N. 1964. Estimate of the genetic load disclosed bv inbreeding. Genetics 
50: 527-529. 

Kosuda, K. 1972. Synergistic effect of inbreeding on viability in Drosophila virilis. 
Genetics 72: 461^168. 



60 



Journal of the Lepidopterists’ Society 



Levene, H., I. M. Lerner, A. Sokoloff, F. K. Ho & I. R. Franklin. 1965. Genetic 
load in Tribolium. Proc. Natl. Acad. Sci. (U.S.A.) 53: 1042-1050. 

Lewontin, R. C. 1974. The Genetic Basis of Evolutionary Change. Columbia Univ. 
Press, New York. 346 pp. 

Malogolowkin-Cohen, Ch., H. Levene, N. P. Dobzhansky & A. Solima-Simmons. 
1964. Inbreeding and the mutational and balanced loads in natural populations of 
Drosophila willistoni. Genetics 50: 1299-1311. 

Morton, N. E., J. F. Crow & H. J. Muller. 1956. An estimate of the mutational 
damage in man from data on consanguineous marriages. Proc. Natl. Acad. Sci. 
(U.S.A.) 42: 855-863. 

Oliver, C. G. 1972a. Genetic and phenotypic differentiation and geographic distance 
in four species of Lepidoptera. Evolution 26: 221-241. 

1972b. Genetic differentiation between English and French populations of the 

satyrid butterfly Pararge megera. Heredity 29: 307-313. 

1977. Genetic incompatibility between populations of the nymphalid butterfly 

Boloria selene from England and the United States. Heredity 39: 279-285. 

1978. Experimental hybridization between the nymphalid butterflies Phy- 

ciodes tharos and P. campestris montana. Evolution 32: 594-601. 

1979. Experimental hybridization between Phyciodes tharos and P. batesii 

(Nymphalidae). J. Lepid. Soc. 33: 6-20. 

Owen, D. B. 1962. Handbook of Statistical Tables. Addison-Wesley, Reading, Mass. 
580 pp. 

PlSANl, J. F. & W. E. Kerr. 1961. Lethal equivalents in domestic animals. Genetics 
46: 773-786. 

Richmond, R. C. 1972. Enzyme variability in the Drosophila willistoni group. III. 
Amounts of variability in the superspecies D. paulistorum. Genetics 70: 87-112. 

Shapiro, A. M. 1978. A new weedy host for the buckeye, Precis coenia (Nymphali- 
dae). J. Lepid. Soc. 32: 224. 

Sittmann, K., H. Abplanalp & R. A. Fraser. 1966. Inbreeding depression in Jap- 
anese quail. Genetics 54: 371-379. 

Sorenson, F. 1969. Embryonic genetic load in coastal Douglas-fir, Pseudotsuga men- 
ziesii var. menziesii. Amer. Natur. 103: 389-398. 

Stone, W. S., F. D. Wilson & V. L. Gerstenberg. 1963. Genetic studies of natural 
populations of Drosophila: Drosophila pseudoobscura , a large dominant popula- 
tion. Genetics 48: 1089-1106. 

Vawter, A. T. & P. F. Brussard. 1975. Genetic stability of populations of Phyciodes 
tharos (Nymphalidae: Melitaeinae). J. Lepid. Soc. 29: 15-23. 



